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Abstract. We present J-H-K' photometry for a sample of 45 high redshift quasars found by the Sloan Digital 
Sky Survey. The sample was originally selected on the basis of optical colors and spans a redshift range from 3.6 
to 5.03. Our photometry reflects the rest-frame SED longward of Lya for all redshifts. The results show that the 
near-IR colors of high redshift quasars are quite uniform. We have modelled the continuum shape of the quasars 
(from just beyond Lya to ~ 4000 A) with a power law of the form i^, oc u°' , and find (a) = —0.57 with a scatter of 
0.33. This value is similar to what is found for lower redshift quasars over the same restframe wavelength range, 
and we conclude that there is hardly any evolution in the continuum properties of optically selected quasars up to 
redshift 5. The spectral indices found by combining near-IR with optical photometry are in general consistent but 
slightly flatter than what is found for the same quasars using the optical spectra and photometry alone, showing 
that the continuum region used to determine the spectral indices can somewhat influence the results. 
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1. Introduction 

The last few years have seen a dramatic increase in the 
number of high redshift optically selected quasars by sev- 
eral surveys (e.g. Warren et al. 1994, Kennefick et al. 
1995). Most prominent amongst these sm'veys is the Sloan 
Digital Sky Survey (SDSS - Fukugita et al. 1996, Gunn 
et al. 1998, Hogg et al. 2001, York et al. 2000, Smith 
et al. 2002, Stoughton et al. 2002, Pier et al. 2002, 
Richards et al. 2002), which has amongst its scientific 
aims, the construction of the largest sample of quasars 
ever, with more than 10^ objects spanning a large range 
of redshift and luminosities. The SDSS has discovered an 
unprecedented number of new high redshift quasars, in- 
cluding more than 200 new quasars at z>4 (e.g. Fan et al. 
2001a, 2000, 1999; Zheng et al. 2000, Anderson et al. 
2001, Schneider et al. 2001) and the most distant quasar 
known to date at z=6.4 (Fan et al. 2003, in press). These 
high redshift quasars have been efficiently selected by their 
distinctive position in color-color diagrams, with charac- 
teristic colors due to the main features of the quasar spec- 
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tra, viz., the power law continuum, the strong Lya emis- 
sion line, the Lya forest absorption and Lyman limit, all 
of which move quasars away from the stellar color locus. 
In this paper we aim to determine the continuum prop- 
erties of high redshift quasars at optical restframe wave- 
lengths. To do this we have obtained J-H-K' photometry of 
a large color-selected sample of high redshift quasars found 
by SDSS. A good knowledge of the quasar continuum 
shape at near-UV/optical restframe wavelengths is impor- 
tant for several reasons. Over the near-UV/optical wave- 
length range the shape of the continuum is usually approx- 
imated with a single power law of the form /{v) oc i^". 
The mean value of a, the dispersion of this value, and 
even the validity of the power law parameterization as 
well as the evolution of such parameters with redshift are 
still under debate (e.g. Vandcn Berk et al. 2001 and refer- 
ences therein). The continuum slopes of quasars are blue, 
with a mean canonical spectral index of a = —0.7 (e.g. 
Richstone and Schmidt 1980), but several results of the 
past few years point to flatter continua (e.g. Francis 1996, 
Natali et al. 1998, Vanden Berk et al. 2001), and indicate 
a different slope at different restframe wavelength ranges. 
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On the other hand the most recent results of Fan et al. 
(2001a) and Schneider et al. (2001), who find steeper av- 
erage indices (respectively a = —0.8 and —0.9) for very 
high redshift quasars, seem to point to an evolution of the 
continuum properties of quasars with redshift. However, 
as discussed in Schneider et al. (2001), for high redshift 
quasars, determination of the continuum slope from opti- 
cal data alone {Xobs < 1^) must rely on a small restframe 
wavelength region with Xresi < 1800 — 1900 A. For ex- 
ample in objects at z~4, the restframe V-band emission 
is shifted to K-band. Therefore near-IR data are essential 
to unambiguously determine the continuum properties of 
high redshift quasars in a way comparable to low-redshift 
objects, so as to get an unbiased measure of evolution. 
A knowledge of the optical continuum shape and its pos- 
sible redshift evolution, is not only important to under- 
stand the quasars itself but also for selecting objects at 
even higher redshifts. Optical colors alone begin to be less 
efficient for selection purposes at z~ 5, since the quasars 
evolutionary track crosses the locus of very low-mass, late- 
type stars. Already, the highest redshift quasars discovered 
to date have been selected by adding J-band photometry 
to the SDSS colors, to distinguish quasar candidates from 
stars (see in particular Zheng et al. 2000 and Fan et al. 
2001a). 

In addition, determining the flux decrement in the Lya 

forest needs an estimate of the contimmm at A < X^ya, 
which is usually an extrapolation from A > 1/^m. If this 
is done by using the classical a = —0.7 index it can lead 
to uncertainties of 5-10% in the computation of the con- 
tinuum decrement. Also for high redshift quasars the ob- 
served optical fluxes must be extrapolated to obtain the 
optical luminosity M^, e.g. to determine the luminosity 
function. Any change of a with redshift could substan- 
tially influence the inferred luminosity. Only with near-IR 
photometry can Mb be determined directly. Lastly, for the 
SDSS quasar search, a knowledge of the continuum slope 
is essential for modeling the SDSS quasar selection func- 
tion which is then used to derive the luminosity function 
(Fan et al. 2001b). 

2. Sample selection 

High-redshift (z> 3.6) quasar candidates were selected us- 
ing color cuts that separate them from the stellar locus 
(e.g. Fan et al. 2000, 2001a). For the color-selected sta- 
tistical sample presented in §4, the color criteria were ap- 
plied after correcting for Galactic extinction using the red- 
dening map of Schlegel, Finkbeiner & Davis (1998). The 
criteria are as follows: 

1. gri candidates, selected principally from the g* — 
r* ,r* — i* diagram: 

(a) i* < 20.05 

(6) u* -g* > 2.00 or u* > 21.00 

(c) g* -r* > 1.00 (1) 

(d) r* -i* < 0.42(g* - r*) - 0.31 or g* - r* > 2.30 

(e) i* -z* < 0.25 



2. riz candidates, selected principally from the r* — 
— z* diagram: 

(a) i* < 20.20 
lb) u* > 22.00 

(c) g* > 22.60 (2) 

(d) r* -i* > 1.00 

(e) i* -z* < 0.47(r* - i*) - 0.48 

The intersections of those color cuts with the g* — 
r*,r* — i* and r* — i* ,i* — z* diagrams are illustrated e.g. 
in Fan et al. (2000) and (2001a). 

Within these color magnitude boundaries, two subsam- 
ples were selected for observations in Spring and Fall. The 
first is the color selected sample in the Fall Equatorial 
Stripe, presented in Fan et al. (2001a), which comprises 
39 objects. The second is a complete sample in the Spring 
Equatorial Stripe, consisting of 55 objects whose redshifts 
have been reported in different papers (Schneider et al. 
2001, Fan et al. 2000, Fan et al. 2003 in preparation , 
Anderson et al. 2001). Since the two samples have been 
selected with the same optical color criteria they can be 
merged together and form a large color selected sample, 
spanning a redshift range from 3.60 to 5.03. 
The completeness of the Fall sample has been exten- 
sively discussed in Fan et al. (2001a) and is around 80 
%, depending slightly on redshift; the completeness for 
the Spring sample should be very similar. As we detail 
below, we have observed a random subsample of 45 of the 
94 quasars. Therefore our sample can be viewed as a sta- 
tistical sample with the above color cuts. 

3. Observations and data reduction 

Observations were carried out in service mode with the 
MAGIC near-IR camera on the 2.2m telescope located at 
Calar Alto (Spain). MAGIC is equipped with a Rockwell 
256x256 pixel NICM0S3 detector array. In the high res- 
olution mode it provides a pixel scale of 0.64" pixel"^ 
and a total field of view of 164" x 164". The observations 
log is shown in Table 1. In this paper we only present 
data from nights with photometric or nearly photomet- 
ric conditions (clear nights, with very stable conditions). 
A fraction of the non-photometric time was used to ob- 
tain J-band snapshots of candidate high redshift quasars 
from the SDSS, selected as i-band dropouts, for which sub- 
sequent spectroscopy was carried out at other telescopes 
(e.g. Fan et al. 2001b). 

Each target was observed through the J, H, and K' 
filters. The K' band (Wainscoat & Cowie 1992) was pre- 
ferred to the K band to reduce the effects of thermal emis- 
sion from the sky and the telescope. Each object was ob- 
served for a total on-target exposure time between 10 and 
12 minutes in each filter. For a satisfactory background 
subtraction, the images were dithered on a 6 position pat- 
tern (4 positions for the first run), with offsets of around 
30" between each frame. We took 6, 12 or 18 frames per 
object depending on the night and on the brightness of 
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Run 


da-t6 


S66ing 


nights /phot 


t 7 


t/f 


tx 


no of obj 


(1) 


(21 


(3) 




(5) 


("1 


(7) 


(8) 


1 


Nov 1999 


1 


1/1 


10m (10sx3) 


10m (10sx3) 


10m (6sx5) 


8 


2 


Aug 2000 


0.9-1.1 


7/3 


12m (10sx6) 


12m (10sx6) 


12m (6sx5) 


19 


3 


Feb 2001 


1< 


4/2 


12m (15sx8) 


12m (15sx8) 


12m (5sxl2) 


17 


4 


Mar 2001 




6/0 


12m (10sx6) 


12m (10sx6) 


12m (10sx6) 





5 


Jun 2001 


1.4-1.5 


2/1 


12m (10sx6) 


12m (10sx6) 


12m (5sxl2) 


6 


6 


Nov 2001 


1-1.5 


3/2 


12m (10sx6) 


12m (10sx6) 


12m (5sxl2) 


12 


7 


Feb 2002 


1.3-1.5 


1/1 


12m (10sx6) 


12m (10sx6) 


12m (5sxl2) 


3 



(1) Run identification (2) Observation date (3) average seeing (FWHM) during the night (4) Total number of nights/ total 
number of photometric-good nights (5) Total integration time in J-band with sub-integration time and number of frames at 
each position (6) The same for H-band (7) The same for K-band (9) Total number of objects observed (including repetition and 
objects that are not in the present sample) 



the object. At each position the total integration time 
was split into short sub-frames to avoid saturation (the 
frames were directly integrated at the telescope) . In Table 
1 we report all these observational parameters, namely the 
total integration time, the sub-integration times and the 
number of frames at each position. 

During each night, 4 or more different standard stars 
(UKIRT faint stars from Elias et al. 1982, complemented 
by those observed in the same field by Hunt et al. 1998) 
spanning a range of colors were observed in each of the 
filters at regular intervals (at least 3 times per night). In 
most cases the standards span the same range of airmasses 
as the targets. From the fall sample a total of 29 objects 
(out of 39) were observed under photometric or nearly 
photometric conditions. From the spring sample a total 
of 16 (out of 55) objects were observed. The selection of 
the objects was only determined by chance and observ- 
ability so it is unbiased with respect to quasar properties. 
Therefore a total of 45 objects form a statistically well 
defined sample with complete ncar-IR photometry. 
Note that the J, H and K' data for a given object were 
taken on a single night, while the near-IR data were ob- 
tained one to two years after the optical fluxes were mea- 
sured from SDSS, so we cannot exclude the possibility 
that variability could influence some of the results, such 
as the determination of the spectral indices of the contin- 
uum slopes. For example, a variability of 10-20% in the 
optical flux could change the measured slope by ~ 0.2. 

Data reduction was performed with IRAF^ . The sky- 
background was subtracted using a composite sky frame 
made for each image from the previous and following im- 
ages. Typically five neighboring frames were used to cre- 
ate a sky image, although in a few cases when the sky was 
changing more rapidly we used only three frames. Dome 
flats were obtained by taking exposures with dome lamps 
on and off and subtracting one from the other; these were 
applied after sky subtraction (an alternative flat field was 
also created by averaging a large number of sky-frames, 

^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the AURA, Inc under coop- 
erative agreement with the NSF. 



masking out the brightest stars). The background sub- 
tracted, flatfielded images were shifted to a common ref- 
erence, with shifts derived from centering the position of 
as many point sources as were visible in the frames. In 
most cases at least four or more such reference stars were 
present, but in a few objects only one reference could be 
used. The registered images were then co-added using the 
IRAF/avsigclip rejection algorithm. With the large num- 
ber of frames for each target (typically 12), it was easy to 
perform cosmic ray rejection. 

Photometric calibration was done using at least five 
standard stars to derive the zero point of the three bands. 
During each good night the scatter in the photometric 
zero point was less than 0.05 (slightly higher for J band), 
without color correction so all the standard stars in a 
given band were simply averaged to give the final cali- 
bration. The photometry is accurate to 0.1 magnitudes 
for the bright part of the sample. Magnitudes were de- 
rived inside circular apertures with radius of 4". The K' 
filter, as opposed to the standard K filter, docs not have 
many published values for standard star photometry, but 
for some stars interpolated K' magnitudes were available 
from the web page developed by Dave Thompson^. For 
those stars which did not have K' magnitudes available 
we simply used the K band magnitude from Hunt et al. 
(1998). 

4. Results and discussion 

4.1. Colors 

The resulting J, H and K' magnitudes with their cor- 
responding errors are presented in Table 2. The photo- 
metric errors are dominated by the sky noise. In Figure 
1 we show the distribution of J, H and K' magnitudes 
with redshift for all quasars in the sample. In Figure 2 
we show their colors (z— J, J— H and J— K') again plotted 
as a function of redshift. Note that the z band values are 
from SDSS photometry, based on the AB magnitude sys- 

^ http://www.astro.caltech.edu/mirror/keck/realpublic/inst/ 
nirc/UKIRTstds.html 
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Fig. 1. The J,H and K magnitude of the quasars plotted 
as a function of redshift, with photometric uncertainties 
indicated by the errorbars. 
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Fig. 2. The near-IR colors of the sample quasars plotted 
as a function of redshift. The dashed lines are the color 
expected from an object with a power law continuum of 
slope a = —0.5 (below) and a = —1.0 (above), while the 
solid lines are the colors obtained by redshifting the com- 
posite SDSS quasar spectrum (Vanden Berk et al. 2001). 

tem, while the IR-colors are on the Johnson (Vega based) 
system. In Figure 3 we present the shape of the spectral 
energy distributions (SED) for all objects, combining the 
SDSS five-band photometry (from Fan et al. 2001a; 2000 
and Richards et al. 2001) with our near-IR colors. The 
observed flux points have been shifted to the restframe 
wavelength and normalized to have the same i-band mag- 
nitude (the i-band flux being of higher S/N than the z 
band). For these plots we have consistently converted the 



near-IR photometry into AB magnitudes using the fol- 
lowing zero points: Jab = J + 0.89, ^ab = H + 1.38 and 
^'ab = ^' + 1-84. These offsets were derived assuming 
that the flux density from Vega is constant within any 
filter and equal to 1620 , 1020 and 688 Jy: the first two 
values are reported on the UKIRT homepage^, while for 
the K' filter we used an extrapolation from a measurement 
available in the MAGIC homepage*. 

In Figure 2 we have also plotted as dashed lines the col- 
ors expected for a quasar if the continuum were a perfect 
power law with a spectral index a — —0.5 and a = —1.0 
respectively. From the continuum only, the colors ex- 
pected would be z-J=1.08, J-H=0.69 and J-K=1.22 for 
a — —0.5. The solid lines in each panel represent the 
colors produced by shifting the composite SDSS quasar 
spectrum by Vanden Berk et al. (2001) to the different 
redshifts. This spectrum is best represented by a slightly 
fiatter power law continuum {a = —0.46) and of course 
contains the contribution of all emission lines and features. 
The two lines deviate substantially from each other only 
in the z— J color for objects at redshift >4.5, due to the 
contribution of CIV entering the z-band, and in the J— H 
color for redshift less than ~ 4, due to the contribution of 
the Mgll and the Fell complex in J band (Richards et al. 
2001, Barkhouseand HaU 2001). 

From the plots it is clear that the individual quasars show 
scatter in their colors around the predicted mean color 
lines. The scatter is far in excess of the measurement er- 
rors. Such scatter is produced by several factors, includ- 
ing the intrinsic differences of the spectral properties of 
quasars, which can have bluer or redder continua and the 
different relative strength of the emission lines in differ- 
ent objects. Note that in Figure 2, the object at z=4.92 
with colors deviating substantially from the expected ones, 
is a BAL quasar SDSSJ 160501.21-011220.6 (Fan et al. 
2000), which shows quite unusual optical colors (see dis- 
cussion in Hall et al. 2002). Apart from this object and 
SDSS J103432. 72-002702.6, which is a mini BAL, there 
are no other BAL quasars in the sample. 
There are only few near-IR measurements available in the 
literature for quasars at a similar redshift (e.g. Rodriguez 
Espinosa et al. 1988, Bechtold et al. 1994, Francis 1996, 
Zheng et al. 2000), all indicating colors similar to what 
we find. 

Note that Zheng et al. (2000) used z-J< 1.5 and J-K 
< 1.8 as additional constraints to select candidate high red- 
shift quasars, whereas Fan et al. (2001b) used only z— J< 
1.5. We see that actually all except 3 of the objects sat- 
isfy the first constraint, and only one does not satisfy the 
second criterion. Indeed all L and M dwarfs stars found 
by combined SDSS plus 2MASS photometry have z— K > 
2 (Finlator et al. 2000, Leggett et al. 2002) so their col- 
ors are quite different from those of quasars. These results 
confirm that the addition of one single near-IR band in- 
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Fig. 3. The spectral energy distribution of all quasars: all 

near-IR magnitudes have been transformed into the AB 
system (see text for details). The observed values have 
been shifted to the restframe of each object and normal- 
ized to have the same i-band magnitude. The lines indicate 
power laws with spectral indices from -0.25 to -1.0 (from 
the lower to the upper line). The object with colors much 
brighter than the other is 2256-1-0047, one of those whose 
continuum is not well represented by a power-law. It could 
be also a variable quasar. 



formation to the SDSS photometry can greatly improve 
the efficiency of finding high redshift quasars. 

4.2. Continuum slopes 

The continuum slopes were computed using all SDSS 
bands that fall fully rcdward of the Lya line (i.e. exclud- 
ing the band that contains Lya) in combination with the 
three ncar-IR bands. As already noted in the previous sec- 
tion, contamination by line emission can be important in 
some bands. The most important lines that are present 
in the quasar spectra at the restframe wavelengths sam- 
pled by our observations are the CIV (Xrest = 1549 A) 
and CIII] {Xrest = 1909 A) lines for the optical photom- 
etry (i and z-band); the Mgll line {Xrest = 2798A) which 
falls in the J-band for objects with z< 4.0 and in the H- 
band for objects with z> 4.36; finally H/3 {Xrest = 4861 
A) which contaminates the K' band for objects with z< 
3.73. Many other lines also fall in these regions but their 
contributions to broad band fluxes are mostly negligible. 
A special case is the very broad feature that appears be- 
tween 2200-4000 A, known as the 3000 A bump, which 
consists of blends of Fell line emission and Balmer con- 
tinuum emission (e.g. Wills et al. 1985, Richards et al. 
2002). No attempt was made to correct for these features 
since their contribution is hard to disentangle. However, 
we have seen from the color distribution in the previous 
section that this contribution cannot be very large, or at 



least it cancels out, otherwise we would observe larger de- 
partures from the values expected from a simple power law 
continuum. Furthermore we are using a quite large wave- 
length baseline (from 7000 to 22000 A) which means that 
these errors would not have a large effect on the derived 
continuum slope. For example an error of 25 per cent in 
the flux of one of the bands would introduce an error of 
< 0.1 in a (in the sense that if the 25% was due to the flux 
contributed by the Fell bump, the true a would be flatter 
than what we measure). 

For a large fraction of the objects in the sample the 

strength of the CIV line is known, having been mea- 
sured from the low resolution discovery spectra (Fan 
et al. 1999; Fan et al. 2001a; Fan et al. 2003 in prepa- 
ration; Anderson et al. 2001; Schneider et al. 2001). The 
strength of the other lines was not known so we used the 
average equivalent width derived from the combined spec- 
tra of SDSS quasars by Vanden Berk et al. (2001). In par- 
ticular we used EW=23.78 A for CIV (when the real value 
was not available), EW=21.19 A for CIII], EW= 34.95 A 
for Mgll and EW=46.21 A for II/3 (these are all restframe 
values). Note that in most cases the contribution of line 
emission flux to the broad band flux is less than 0.1 magni- 
tudes, so it is comparable to or less than the photometric 
error. 

We subtracted the line contribution from the relevant 
broad band values and then converted magnitudes into 
flux. The corrected continuum fluxes were then modeled 
with a power law of the form oc v" and the best fitting 
spectral index a was derived by least squares minimiza- 
tion. The fit by a power law was considered acceptable 
when the probability that a value of chi-square as poor as 
the value found should occur by chance was larger than 
0.05. For 5 objects in the sample (0019-0040, 0035+0040, 
1625-0001, 2256+0047 and 2306+0108) this requirement 
was not satisfied, so we report the values we have obtained 
in Table 2, but we do not include them in the discussion. 

In Figure 4 we present the distribution of a from the 
combined Fall and Spring sample. The distribution is es- 
sentially confined between 0.2 and -1.1 (with one object at 
-1.8), with an average spectral index of a = —0.57 and la 
dispersion of 0.33. The mean was weighted by the errors 
on the individual a measurements. The median value is 
a = —0.51. The distribution is not perfectly symmetric, 
but slightly skewed towards steeper indices. In Figure 4 
(right panel) we also show for comparison the distribu- 
tion of spectral indices derived by Ivezic et al. (2002) for 
about 6800 SDSS quasars with magnitude brighter than 
i=19, spanning a large range of redshifts (see also dis- 
cussion in Richards et al. 2002). The two distributions 
appear similar. The values found for our sample are con- 
sistent with the average power law index derived from 
the composite SDSS quasar spectra, spanning a redshift 
range 0.04< z < 4.79 (Vanden Berk et al. 2001) which is 
a = —0.47. This index is valid for a large spectral range, 
1216A < A < 5000 A restframe: at longer wavelengths, 
the slope changes to a =-1.58. A similar result was ob- 
tained by Francis (1996) who derived a using photometric 



6 Laura Pentericci et al.: Near-IR properties of high-z QSOs 

Table 2. Near-IR multiband photometry of high redshift SDSS quasars 
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References for redshift: F99 Fan et al. 1999; FOl Fan et al. 2001a; F03 Fan et al. 2003 in preparation; AOl Anderson et al. 
2001; SOI Schneider et al. 2001. 



estimates, including near-IR photometry, for a sample of 
quasars spanning a redshift range 0.33 to 3.67 and with 
a median redshift of 2. In Figure 4 wc plot the spectral 
indices derived by Francis (1996) for all the objects with 
redshift above 2: for these objects the average spectral in- 
dex is a = -0.49 (median -0.43). Natah ct al. (1998) find 
a similar result for a sample of bright quasars with red- 
shift up to 2.5, with a=-0.65. However this slope is based 
on the continuum shortward of the 3000 A bump, with a 



much flatter index at 4000 A. One discrepant result comes 
from Kuhn et al. (2001), who derived the slopes of a sam- 
ple of bright quasars with redshift around 3, using both 
photometry and spectroscopy. They find slightly flatter 
average spectra, with a = —0.34 (median -0.29); however 
wc note that their spectral slopes for the z^O comparison 
sample of quasars are also flatter than what was found in 
other studies. 
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a 

Spectral Index 



Fig. 4. Left: the distribution of the quasars continuum power-law index a from the combined fall and spring samples. 
The dotted line represents the distribution of spectral indices from a sample of quasars at redshift between 2 and 3.5 
from Francis (1996). Right: the distribution of spectral indices derived by Ivezic et al. (2002) for about 6800 SDSS 
quasars with magnitude brighter than i=19. The a distribution for a subsample of 440 radio-loud quasars, indicated 
by the triangles, is skewed towards more negative values. 



In Figure 5 we plot the derived values for the Fall sam- 
ple, and compare them to the spectral indices derived for 
the same quasars in Fan et al. (2001a) from the low res- 
olution spectra and optical photometry. The error bars 
plotted in the figure indicate the Itr error. Basically all 
spectral indices are consistent with those derived from 
optical data within 3a. However the indices derived from 
the combined near-IR and optical broad band photome- 
try tend to be flatter (bluer) than those derived from the 
optical data alone and the dispersion of a values around 
the mean is also somewhat smaller. The average spectral 
index derived from the photometry for the 29 objects of 
the fall sample alone is a = —0.59 ± 0.34 with a median 
value of —0.47, whereas for the compete fall sample (all 
39 quasars) Fan et al. (2001a) derived an average value of 
a = —0.88 ±0.36 (note however that they apphed weights 
that included the error measurement as well as the de- 
tection probability function). Similarly Schneider et al. 
(2001) derived an average a = —0.93 ± 0.31 from spectra 
of high redshift SDSS quasars, some of which are included 
in our sample. 

The apparent discrepacy probably arises because the in- 
dices are derived from different wavelength regions of the 
continuum emission: the indices derived by Fan et al. 
2001a and by Schneider et al. (2001) are from the con- 
tinuum immediately blueward of the Lya up to 9000 A 
(or less) observed emission, and thus include only a few 
hundred A in the restframe (1250 to 1600-1800 A depend- 
ing on the redshift). Furthermore, part of the Fe complex 
in the 1500 - 2000 A region could make the i and z band 
brighter, thus producing a steepening of the spectral in- 



dices. On the other hand our photometry spans a much 
larger wavelength region and, most important, includes 
longer wavelengths: at redshift 3.65 (our lower redshift) 
the photometry cover from 1500 to 4700 A (restfr ame 
whereas at redshift 5, our most distant object, the cover- 
age goes from ~ 1150 A up to ~ 3700 A. Fig 3 illustrates 
that the continuum immediately longward of Lya is some- 
what redder than at longer wavelengths. 

We conclude that the continuum properties of the high 
redshift quasars in our sample at optical restframe wave- 
lengths are comparable to those of their lower redshift 
counterparts, with no significant change with epoch. We 
also find that the spectral slopes change somewhat de- 
pending on the wavelength regions used to measure them, 
as already indicated by the results of Vanden Berk et al. 
(2001) and Natah et al. (1998). 

5. Summary 

We have presented near-IR photomety (J,II,K') of a sam- 
ple of 45 color-selected SDSS quasars at 3.6< z <5. We 
have determined the slopes of the optical continuum at 
restframe wavelengths 1200 A< \rest < 5000 A, and we 
find an average spectral slope a = —0.57 and a 1 tr dis- 
persion of 0.33. This value is consistent with those found 
by many authors for lower redshft quasars at similar rest- 
frame wavelengths, so we conclude that there is no ev- 
idence for a redshift evolution of the quasar continuum 
properties, as it was previously suggested. By comparing 
our slopes with those determined by Fan et al. (2001) for 
the same objects, using the optical spectroscopy and pho- 
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Fig. 5. The spectral index of the continuum emission: on 
the y axis is the value derived by Fan et al. 2001 from the 
low resolution spectra, on the x-axis the value derived from 
optical and near-IR photometry. The error bars represent 
the 1 fj uncertainties. 

tometry alone, we can see that the addition of near-IR 
data yields a somewhat bluer continuum and more ro- 
bust evaluation. This implyes that the determination of 
the spectral slope changes somewhat depending on the 
wavelength used to measure them, as already indicated 
by previous results. 
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